Nonextensive hydrodynamics for relativistic heavy-ion collisions 
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The nonextensive one-dimensional version of a hydrodynamical model for multiparticle produc- 
tion processes is proposed and discussed. It is based on nonextensive statistics assumed in the form 
proposed by Tsallis and characterized by a nonextensivity parameter q. In this formulation the 
parameter q describes some specific form of local equilibrium which is characteristic for the nonex- 
tensive thermodynamics and which replaces the local thermal equilibrium assumption of the usual 
hydrodynamical models. We argue that there is correspondence between the perfect nonextensive 
hydrodynamics and the usual dissipative hydrodynamics. It leads to simple expression for dissipa- 
tive entropy current and allows for predictions for the ratio of bulk and shear viscosities to entropy 
density, £/s and n/s, to be made. 

PACS numbers: 24.10.Nz, 25.75.-q 
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I. INTRODUCTION 

Multiparticle production experiments are our main 
source of information on multiparticle production pro- 
cesses in which the initial kinetic energy of two projectiles 
is to a large extent converted into a multitude of observed 
secondaries. This is especially true in the case of heavy 
ion collisions in which one expects the formation of a 
new hadronic state of matter, the Quark Gluon Plasma 
(QGP) [1]. Such processes call for some form of statisti- 
cal approach, which is usually based on the Boltzmann- 
Gibbs (BG) statistics. On the other hand, in the case 
of a multiparticle production process, conditions lead- 
ing to BG statistics are satisfied only approximately at 
best. This is because, among other things, hadronizing 
systems experience strong intrinsic fluctuations and long 
range correlations 0, H[ , which can be interpreted as sig- 
nals of some dynamical, nonequilibrium effects showing 
up [J]. It is therefore difficult to expect the occurrence 
of the usual (local) thermal equilibrium, instead one has 
some kind of stationary state. It turns out that these 
phenomena can be incorporated, at least to some extent 
and without going into deeper dynamical considerations 
concerning the sources of such fluctuations, in the formal- 
ism of the nonextensive statistics (which we shall apply 
here in the manner proposed by Tsallis, see ||) in the 
form of a more general equilibrium summarily described 
by a single parameter q [s 0, fl ■ This parameter char- 
acterizes the corresponding Tsallis entropy, S q , which in 
such approach replaces the usual BG entropy to which 
it converges for q — ► 1. Because such systems are in 
general nonextensive, the parameter q is usually called a 
nonextensivity parameter. 

Such approach has been successfully applied to mul- 



tiparticle production processes, both by using nonex- 
tensive versions of the respective distribution functions 
PS El El El, EI EE d, or by deriving such distribu- 
tions from the appropriately modified nonextensive ver- 
sion of the Boltzmann transport equation [1, 0, H E3j EH • 
In both cases, this amounts to replacing the usual expo- 
nential factors by their ^-exponential equivalents, 

P BG (E) =exp(-£/T) P q (E) = ex Pq (-E/T) 

= [l-(l-q)E/T] 1 ^. (1) 

Notice that Pbg(E) = P q =\(E). In all these applications 
one finds that q > 1 . It represents effect of some intrinsic 
fluctuations existing in the hadronizing system and re- 
vealing themselves as fluctuations of its temperature or 
of the mean multiplicity of secondaries [57| . Generically 
one has that [2j 
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where, depending on the kind of process considered, T 
can be replaced by some other variable [H EH EH . Be- 
cause different observables are sensitive to different kinds 
of fluctuations, it is natural to expect that the y a re de- 
scribed by different values of the parameter q [14| . For 
example, single particle, one-dimensional distributions in 
longitudinal phase-space, like dN/dy, are most sensitive 
to fluctuations of the mean multiplicity (n) of secondaries 
(l3| and are described by the nonextensivity parameter 
q = qi, of the order of q^ — 1 ~ 0.1 -j- 0.2 [58]. However, 
distributions in transverse momenta, dN/dpr, which are 
believed to probe the local thermal equilibrium of the 
hadronizing system (assuming that such a phenomenon 
indeed occurs) and serve as a source of information on 
the temperature T of the hadronizing system, are very 
sensitive to fluctuations of this temperature represented 
by the nonextensivity parameter q = qT [HI ]" changes as 
small as q^ — 1 ~ 0.01 -j- 0.05 already substantially affect 
the resultant pt spectra [5j| Notice that data indicate 
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that ql — 1 3> 9t — 1) i- e -i fluctuations governed by 
arc dominant and therefore for the whole system q ~ 
0. 

Among statistical approaches to the multiparticle pro- 
duction processes a specially important role is played by 
hydrodynamical models [21| . which so far are all based 
on the BG statistics. The existing general nonexten- 
sive version of fluid dynamics discussed in |22j is not 
suitable for applications to multiparticle production pro- 
cesses (among other things because of its non covari- 
ant formulation). We would like to fill this gap and 
present a fully covariant hydrodynamical model based 
on q-statistics which can be applied to a multiparticle 
production processes, especially to relativistic heavy- ion 
collisions. Because of exploratory character of our pa- 
per we limit ourselves to the (1 + 1) dimensional case 
only and confront our results with rapidity and trans- 
verse momenta distributions obtained recently at RHIC 
leaving the most detailed studies of all aspects of avail- 
able experimental data for future investigations. 

The hydrodynamical model of multiparticle produc- 
tion means, in fact, a number of separate problems con- 
nected with the consecutive steps of the collision pro- 
cess: the choice of initial conditions summarizing the 
preparatory stage of collision (it should end in some form 
of local thermal equilibrium), the choice of equation of 
state (EoS) of the quark-gluon and/or hadronic matter 
being equilibrated, further hydrodynamical evolution of 
this matter assumed to form a kind of fluid, final con- 
version of this fluid into observed secondaries. Because 
dynamical factors underlying each step are different the 
resulting fluctuation patterns can also differ, presumably 
leading to the parameter q changing during the collision 
process. However, in the present study we shall restrict 
ourself only to the case of nonextensive parameter q re- 
maining the same in the whole collision process [6f3 |. 

Recently there is renewed interest in dissipative hydro- 
dynamical models [H [M [H, [H, HE HI Hi, M, M, 11 , 
which is prompted by the apparent success of hydro- 
dynamical models in describing data obtained at RHIC 
[31 BUI and by the recent calculations of transport co- 
efficients of strongly-interacting quark-gluon system us- 
ing the AdS/CFT correspondence [36[. The question 
addressed is whether, and under which circumstances, 
dissipative hydrodynamics is really needed and how it 



should be applied. The reason is that formulation of 
the relativistic hydrodynamic equations for dissipative 
fluid, suffer from ambiguities in the form they are writ- 
ten [27] . unphysical instability of the equilibrium state 
in the first order theory [28[ and the loss of causality 
in the first order equation approach [3l|, to mention a 
few. In this work we argue that there is a link be- 
tween dissipative hydrodynamics (^-hydrodynamics) and 
nonextensive hydrodynamics (q-hydrodynamics) we are 
proposing, which we call nonextensive/ 'dissipative cor- 
respondence (NexDC). In particular, in Section [V] we 
demonstrate that it is possible to write some of the cor- 
responding transport coefficients of the produced matter 
(believed to be quark-gluon plasma (QGP)) as (implicit) 
functions of the nonextensivity parameter q. The merit of 
using the NexDC is that one can formulate and solve the 
q- hydrodynamic equations of perfect nonextensive hydro- 
dynamics (or perfect q-hydrodynamics) in analogous way 
as for the usual perfect hydrodynamics, which seems to 
be a priori much easier task. Although this does not 
fully solve the problems of d- hydrodynamics, neverthe- 
less it allows us to extend the usual perfect fluid approach 
(using only one new parameter q) well behind its usual 
limits, namely toward the regions reserved for dissipative 
approach only. 

The paper is organized as follows. We start in Section 
im with a short reminder of the nonextensive version of 
kinetic theory from which nonextensive hydrodynamics 
is derived in Section Hill Section HVl contains examples of 
comparisons with experimental data, whereas in Section 
IVl we discuss the possible physical meaning of the pro- 
posed q-hydrodynamics. We end with Section IVII which 
contains our conclusions and summary. Some specialized 
topics and derivations are presented in Appendices [XI [Bl 
[0 and El 



II. RELATIVISTIC NONEXTENSIVE KINETIC 
THEORY 

Following [ItJ we start with nonextensive version 
of Boltzmann equation (the metric used is: g^ v — 
diag(l, -1,-1,-1)), 



C q (x,p)-- 



C q (x,p) 

1 rd 3 Pl d 3 P ' d 3 P [ 



pS p'° Pi 



^{M/^i]wV,pilp,Pi) - M/^iMp-piIp'^'i)}- 



(3a) 
(3b) 



Here f q (x,p) and C q (x,p) are q- versions of the, respec- transition rate between the two particle state with ini- 
tively, corresponding phase space distribution function tial four-momenta p and pi and some final state with 
and the q-collision term in which W(p,pi\p' ,p'x) is the four-momenta p' and p[ whereas h q [f q , f q {\ is the corre- 
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lation function related to the presence of two particles in 
the same space-time position x but with different four- 
momenta p and pi , respectively. Notice two distinct fea- 
tures of Eq. (|3"a|) : (i) it applies to /| (= (f q ) q ) rather 
than to f q itself and (ii) in C q one assumes a new, q- 
generalized, version of the Boltzmann molecular chaos 
hypothesis @, 0, i, S H3 according to which 



Kifqi fqi] = ex P 9 t ln 9 fq + ln ? -M 



(4) 



where exp q (X) = [1 + (1 - and \n q (X) = 

[X^-d - l] /(l - <?)). Eq. Q is our central point, it 
amounts to assuming that, instead of the strict (local) 
equilibrium, a kind of stationary state is being formed, 
which also includes some interactions (see [(1 0, [M Q ) ■ 

With such correlation function one finds that diver- 
gence of the entropy current, which we define as 



< y^^0"{/g(>,P) fo q f q (x,p) - f q {x,p)}, 



(2 



(5) 



r 



is always positive at any space-time point, 
d^(x) > 0, 



(6) 



(this fact is equivalent to demonstrating the validity 
of the relativistic local if -theorem when using this q- 
entropy current). 

To get explicit form of the distribution functions 
f q (x,p) we proceed now as follows 17J. At first, using 
momentum conservation condition in two particle colli- 



sions, p^ 
invariant: 



Pi 



Pi , we form following collision 



^TT 1p(x,p) C q (x,p) 
p" 



0. 



(7) 



where ip(x.p) — a(x) + b fJ- (x)p tJ ' with arbitrary func- 
tions a(x) and b^(x). We assume here that the 
correlation function h q is symmetric and positive, 
h q [f,fi] = hq[fi, f] > 0, and that detailed balance 
holds, W(p,Pib',K) = W(p',p[\p, Pl ). For a(x) = 
and b^{x) = constant one gets the q- version of the local 
energy-momentum conservation [17] . 



d v Tf{x) = 0, 



(8) 



with a nonextensive energy-momentum tensor defined by 

fd 3 p 



0) one 



At the same time for a(x) = constant and b ll (x 

gets m 



d 3 p p 



(27Tft)3 p0/?( X 'P) 



= 0, 



(10) 



this implies that (dfl stands for the corresponding phase 
space volume element) 



d_ 

dt 



Jdn f q (x, P ) = o } 



(11) 



i.e., that the normalization Z q = Jdfl f q {x,p) is con- 
served as well [6l| . Since the divergence of the q-entropy 
current can also be expressed via collision invariant, 



f) s m — \ 

_ (27Tft)3 



F[\n q f q (x,p)}, 



(12) 



demanding that d^s^x) = one finally obtains 



f q (x,p)= l + (l-q)(a(x) + b li (x)p») 



1/(1-9) 



> (13) 



which represents the distribution function in a stationary 
state. Setting a(x) = and ^(x) = —uji(x)/knT q (x) 
(where T q (x) is the temperature function [62] ) one ob- 
tains the well known (unnormalized) Tsallis distribution 
function, 



fq{x,p) = 



l-(l-9) 



V^ q 



1 1/(1-9) 



kBT q (x) 



exp 



k B T q (x) 



(14) 



where ti^ (a:) should be regarded as a hydrodynamical flow 
four vector (hereafter we use the convention that h — 
k B = c = 1). 

We shall now assume that the g-modified energy- 
momentum tensor T£ v can be decomposed in the usual 
way in terms of the g-modified energy density and pres- 
sure, e q and P q , by using the g-modified flow (such 
that for q — ► 1 it becomes the usual hydrodynamical flow 
w M and in the rest frame of the fluid u{f = (1, 0, 0, 0)), 



Tf» = ( £g + p g )«x- V 

— tgUq a q 1 q^q ■ 



(15a) 
(15b) 



where 



q <r 



Denoting e = p /T, z = m/T, 
with g being degeneracy factor depending on the type 
of particles composing our fluid, one gets that in its rest 
frame (or in q-equilibrium) 
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Sq— Uq^lq Uq^ 


gT q 4 

2tt 2 j 


fj 1 2 

/ de v e £ — 


9 2 n /i 
2 2 e [1 — (1 - 


\ 19/(1-9) 


(16a) 


P = — -TV A 

3 q qfJ,L/ 


gT q 4 

' 2tt 2 


J de\/ e 2 - 


- z 2 e[l - (1- 




(16b) 


Sq — SqUq^ 


git 

2n 2 , 


j de\/ e 2 — 


z 2 e|e[l - (1 


- q)e] q ^ 1 -^ + [1 - (1 - q)e] 1 ^ 1 - q ^\ 


(16c) 



(notice that for q < 1 the integration range is limited to 
z < e < 1/(1 — q) in order to keep the integrand positive). 
It is straightforward to check that, in the baryon free case 
to which we shall limit ourselves here, 



and 



TqSq — £q Pq 



dPq _ 

dT q ~ Sq -> 



(17a) 



(17b) 



i.e., that the usual thermodynamic relations also holds 
for the g-modified quantities. 



III. THE NONEXTENSIVE 
HYDRO DYNAMICAL MODEL - 
(/-HYDRODYNAMICS 

A. Equations of nonextensive flow - q-flow 

Our starting point in formulating the q- 
hydrodynamical model is Eq. ||5J) with energy-momentum 
tensor T^ v given by Eq. Because, due to the q- 

version of thermodynamical relations, Eq. (| 1 T[> . in our 
case Eq. © also implies conservation of g-entropy, 



8 s M 

u H t >q 



0, 



(18) 



with s^(x) defined by Eq. ([5]), which can also be written 
as 



(19) 



we have only one general equation which, when written 
using general coordinates and covariant derivatives [63|, 
takes the form: 



TZ = [(e q + P q )K< ~ .„ = 0. 



(20) 



This means that we are dealing here with perfect q- 
hydrodynamics. 

Before proceeding further some specific points of q- 
hydrodynamics, not mentioned in the general derivation 
oresented in Section [TIJ must be kept in mind. At first 
notice that, whereas in the usual perfect hydrodynam- 
ics (based on BG statistics) entropy is conserved in hy- 
drodynamical evolution both locally and globally, in the 



nonextensive approach it is conserve only locally, cf., Eq. 
(|18[) . The total entropy of the whole expanding system 
is not conserved, because for two volumes one finds 
that 



(21) 



where Sq are the corresponding total entropies. Al- 
though, strictly speaking, the hydrodynamical model 
does not require global entropy conservation but only its 
local conservation, the above feature of g-hydrodynamics 
should be always remembered (the consequences of this 
fact will be discussed in more detail in Section |VJ). Sec- 
ond point concerns the causality problem. To guar- 
antee that hydrodynamics makes sense, there should 
exists some spacial scale L such that volume L 3 con- 
tains enough particles composing our fluid. However, in 
case when there are some fluctuations and/or correla- 
tions with some typical correlation length I, for which 
we expect that I > L, one has to use nonextensive en- 
tropy S q L ' (cf. Eq. I|2ip) and its defined locally density 

s q (x) — S q L '/L 3 . When formulating the corresponding 
g-hydrodynamics one takes, as usual, limit L — > 0, in 
which case explicit dependence on the scale L vanishes 
whereas the correlation length leaves its imprint as pa- 
rameter q. In this sense perfect g-hydrodynamics can be 
considered as preserving causality and nonextensivity q 
is then related with the correlation length I. One can 
argue that very roughly that q ~ l/L e g > 1, where L c g 
is some effective spacial scale of the g-hydrodynamics. 
Note here that if the correlation length I is compatible 
with the scale L e ff, i-e., I ~ L e g, one recovers condition 
of the usual local thermal equilibrium and in this case 
the g-hydrodynamics reduces to the usual (BG) hydro- 
dynamics. 

Let us now continue our presentation. When con- 
tracted with the velocity u qv or with the projection ten- 
sor AqXv = g\ v — u q \u qu it leads to the following two 
equations: 



U^df.Eq + (e, + P q )u 

( £ q + Pg^AgXuU^ 



PqUq v9 % = 0, 



(22) 



-AqXud V Pq - PqA qX „g% = 0. (23) 

These are the equations to be solved now for the (1 + 1) 
dimensional case. We shall assume longitudinal expan- 
sion only and introduce proper time r and the space-time 
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rapidity 77: 



t= yjt 2 - z 2 , 

1, t + z 
77 = — In . 

' 2 f - z 



(24a) 
(24b) 



The corresponding metric tensor in this (r-rf) space is 
guv _ diag^l,— . The corresponding four velocity 
of our fluid can be expressed by the local fluid rapidity 
a a (x) as 



In this case Eq. 



1 



(25) 



d£q + Vq_dSq_ + 

dr t drj 



cosh (a q — 77) , — sinh (a„ — 77) 

T 

reduces to (here v q = tanh(a g — 77)) 

(26) 



whereas Eq. (|23|) reduces to the g-generalized relativistic 
Euler equation (cf. Appendix [Al for details): 



To solve these equations one needs additional input in 
terms of EoS, P q = P q (e q ), and the choice of boundary 
conditions, which we set as v q — at 77 = (because 
of the symmetry a = 0). At 77 = Eqs. (f2"2"| and 
reduce to 



de q 
dr 

dj\ 
drj 



P q da 



drj 



r/=0 



= 0. 



»7=0 



(28a) 
(28b) 



B. Nonextensive Equation of State - g-EoS 

The next important ingredient of any hydrodynamical 
model is an Equation of State (EoS) defining a relation 
between the pressure and the energy density, which de- 
pends on the properties of the hadronic matter under 
consideration. In this work we shall only work with EoS 
for the relativistic pion gas (with m T = 0.14 GeV) with- 
out considering different phases of hadronic matter as 
in 39]. The pressure P q and the energy density e q can 



be connected in the form of EoS, P q (e q ), using Eqs. (|16a[ ) 
and (|16bp . However, differently than in the usual cases of 
9=1, the additional freedom represented by the nonex- 
tensivity parameter q makes P q = P q {e q ) ambiguous and 
one has to additionally specify the possible variations of 
the parameter q during the evolution process. In what 
follows, we shall assume that the parameter q remains 
fixed during the the whole evolution of our hadronic fluid. 
We therefore get different EoS for different (but fixed) 
values of the parameter q, examples of which are shown 
in Figure [TJ It displays the ratio P q /s q as a function of 
energy density e q for different values of q = 1.0, 1.1 and 



0.35 



0.30 



0: 

"s- 0.25 



0.20 




10 10 



s q [GeV/fm 3 ] 



FIG. 1: Illustration of EoS for the relativistic nonextensive 
pionic gas (m = 0.14 GeV): the ratio P q /e q is plotted as 
function of energy density e q for different nonextensivity pa- 
rameters q. Notice that the q dependence of EoS shows mainly 
at low and very low energy densities. 



1.2; the temperature T q was varied in the range 0.1 -7- 500 
MeV. It turns out that the g-dependence is confined only 
to the very low energy density region (supporting there- 
fore the previous results on this matter in [lTlEO])- In the 
region of interest, e q ~ 0.1 -j- 5.0 GcV/fm 3 , the changes 
are very small and rapidly vanish with increasing e. 



C. Nonextensive initial conditions - g-initial 
conditions 

To solve equations of (1 + 1) ^-hydrodynamics one has 
to decide on the initial conditions from which the hy- 
drodynamical expansion starts. They must contain some 
form of the local thermal equilibrium, which - as we as- 
sume - is established during the collision process. Accord- 
ing to recent estimations this can happen very rapidly, 
already in the first 1 fm of expansion, if caused by some 
violent, nonperturbative mechanisms operating at this 
stage j64[. It is thus natural to expect that there must 
also exist some intrinsic fluctuations already present in 
this preparatory stage of collision process which, accord- 
ing to our philosophy, should be accounted for by the 
same ^-statistical approach as that used to form the q- 
hydrodynamics. Following [3!||43| we shall use Gaussian 
initial conditions interpolating between two extreme sit- 
uations, the one described by the Bjorken scaling type 
model |44| and the other corresponding to the Landau 
Model [45], but we shall modify them accordingly by 
changing exp(AT) to exp q (X) (as in Eq. Q, it reduces 
to the usual Gaussian of [39| for q = 1). As in [39| ini- 
tial conditions are imposed for the energy density e q ex- 
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pressed as a function of rapidity 77: 



2ct 2 



(29) 



In what follows we shall also require that the q-fluid and 
the space-time rapidities coincide at tq, 



1(7-0,57) = V- 



(30) 



In all calculations presented in this paper we shall as- 
sume for simplicity that e q and a q are independent of 
the transverse coordinate. However, the remaining two 
parameters, e^ 111 -* and a are not independent because one 
has to reproduce the total energy Etot allocated to the 
fluid which is fixed by the conditions of the experiment, 



E tot = ttA t t / di] e g (r ,?y) 



(31) 



where Aj- is the transverse size of the fluid, r is the 
initial proper time r when the fluid starts to expand. 
The E to t can be obtained knowing the mean number of 
participating nucleons N pa I t a nd the total energy loss per 
participating nucleon AE [65] 



E to t — N paxt AE. 



(32) 



The possible initial conditions vary therefore between 
two extremal situations (cf., Fig. [2]): 



(i) The width a in Eq. (I29| is assumed to be fixed 
and kept constant, but its distribution to vary by 
changing e' ln ) to reproduce the fixed total energy 
E to t when the nonextensive parameter q changes. 

(ii) The maximum energy density e^ m ' in Eq , (|29p is 
assumed to be fixed and kept constant, but its dis- 
tribution to vary by changing a to reproduce the 
fixed total energy E to t when the nonextensive pa- 
rameter q changes. 

We would like to stress at this point that such q- 
dependent initial conditions introduce a completely new 
element to hydrodynamical models, not discussed previ- 
ously. The real situation will interpolate in an a priori 
unknown manner between these two extremes, therefore, 
in what follows, we shall restrict ourself only to them. 
As one can see in Fig. [51 whereas the first extreme in- 
troduces sizable g-dependence, the second one leads to 
only minor effects. In both cases, increasing the value of 
q results, as expected (13|, in the enhancement of tails 
for large values of 77. Following [39j our calculations were 
performed for Au+Au collisions at ^/snn = 200 GeV, 
using results reported by the BRAHMS experiment [46| . 
with Etot = 26.1 TeV , V part = 357, AE = 73 ± 6 GeV 
and with At = 6.5 fm, tq = 1.0 fm, see Table]!] In 
Fig. [3] are shown, the initial conditions for the, respec- 
tively, energy density £5(70,77), entropy density s q (T ,r]) 
and temperature T g (70,17) m the case of a = 1.25 fixed 



Initial condition: n = 1.25 
— ¥ =1.00, £^"'=28.7 GeV/fm 3 
■ - ? =1.05, s£ n) =25.9 GeV/fm 3 
y =1.10, e£ ii) =22.4 GeV/fm 3 
? =1.15, £<| n) =16.8 GeV/fm 3 




Initial condition: E (in) =28.7 [GeV/fm 3 ] 



q =1.00, (T = 1.25 
q =1.05, tr = 1.20 
¥ =1.10, a = 1.15 
q =1.15, <r= 1.08 




FIG. 2: Examples of two different types of initial conditions: 
(a) - type (i) initial conditions with fixed a = 1.25 and e*- 111 ' 
varying according to Eq. (|31|l : (b) - type (ii) initial condi- 
tions with fixed e (in) = 28.7 [GeV/fm 3 ] (this is the value 
corresponding to the q = 1 situation in (a)) and a varying 
according to Eq. (J3TJ). 



(d)-(f)) and reproducing the initial energy E to t = 26.1 
TeV. We start with e q (panels (a) and (d) of Fig. |3J 
which is given by Eq. (12"9")) . use it to solve Eq. (|16a|) and 
find T q (ro) (panels (b) and (e) of Fig. [3]) and eventually 
obtain s q (to , 77) using these results and Eq. (|16c|l (panels 
(c) and (f) of Fig. |3J) . 



D. Examples of q- hydrodynamical evolution of 
different thermodynamical quantities 



Let us now demonstrate examples of g-hydrodynamical 
evolution of different thermodynamical quantities and of 
the fluid rapidity. Calculations were performed using the 
method presented in Appendix [E] In Fig [4] we present 
the evolution of the, respectively, energy density e q , tern- 



(panels (a)-(c)) and — 27.8 GeV/fm 3 fixed (panels perature T q and entropy density s q using initial condi 
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TABLE I: List of parameters of the initial conditions used in Fig. [2] The initial temperature Tin = T(to,tj—0) is shown for two 
types of EoS: for the relativistic nonextensive pion gas for some selected values of q > 1 and for the usual BG pion gas with 
q=l. 



Initial condition: 


cr=1.25 fixed 
















Etot 


£ (in) 


a 


EoS 
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9=1.05 


g=1.10 


9=1.15 


9=1.00 
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FIG. 3: Dependence of the initial conditions on the parameter q. The initial conditions for the energy density e q (ro,ri) (panels 
(a) and (d)), the corresponding temperature T q (ro, 77) (panels (b) and (e)) and the entropy density s 9 (to, rj) (panels (c) and (f)) 
are plotted for different values of q and for two types of initial conditions: (i) when the initial width a = 1.25 remains fixed 
(panels (a)-(c)) and (ii) when the initial maximal energy density e^ ia ' = 27.8 GeV/fm 3 remains fixed (panels (d)-(f)). In both 
cases we assume that E to t is the same for all 9's and equal to E to t = 26.1 TeV whereas to in Eq . (I29|l is put equal to to = 1.0 
fm. Temperatures T q for different values of 9 displayed on panels (b) and (e) are determined by solving Eq. (|16a[) whereas the 
entropy densities s q displayed on panels (c) and (f) are obtained from Eq. (|16c[) using values of T q displayed in panels (b) and 
(e). 



tions discussed in Section IIII CI (the exact values of rel- 
evant parameters for both types of initial conditions are 
listed in Table [11} • One can see that the initial func- 
tional forms of e q (r, 77), T q (r, 77) and s 9 (t, 77) generally 
follow their original Gaussian shapes assumed at the ini- 
tial time To for q = 1.0, 1.05 and 1.1. On the other hand, 
during the whole hydrodynamical evolution both the en- 
ergy density e q and the temperature T q calculated for 
q > 1 are smaller than those for q = 1 for t > tq = 1 
fm (see Table UT|) . That is even true for the initial con- 
dition with fixed a = 1.25, for which the initial energy 
density £9=1.1(70,77 = 3) > e g= i.o (r , r? = 3), in which 
case, after the q-hydrodynamical evolution is completed, 
one observes that e q= \ > e q> \. The same trend is also 
observed for the temperature, i.e., T q= \ > T q> \ for all 



t and 77. However, the corresponding entropy density 
s q evolves differently: for both type of initial conditions 
and any 77, inequality relations between s 9= i.i(t, 77) and 
Sg=i.o( T , v) given at initial r = tq are preserved during 
hydrodynamical evolution. 

In what concerns the fluid rapidity a q , it is always 
set to be equal to a g (T ,77) = 77 at r = To. However, 
the pressure gradient, which is characteristic to Gaussian 
type initial conditions applied here, accelerates the fluid, 
therefore a q evolves with time r (actually, this is true 
even for <7=1), see Fig. [5] In these figures the fluid rapid- 
ity a q (actually its deviation from the rapidity 77, a q — 77) 
is shown as a function of r and the corresponding energy 
density e q . Notice that a q — rj = at To for the whole 77 
space (i.e. for all region of the e q ). As shown in Fig. O 
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TABLE II: The values of energy density e q , entropy density s q and temperature T q at r\ — 0.0 and 77 = 3.0 and at r = 5.0 and 
25.0 fm for q = 1.0 and q — 1.1 for two extremal types of initial conditions: with fixed a = 1.25 (upper panel) and with fixed 
£™=27.8GeV/fm 3 (lower panel). 



Initial condition: (7=1.25 fixed 




e q [GeV/fm ;i ] 


T q [GeV] 


s q [1/fm*] 


r 


1 fm 


5 fm 


25 fm 


1 fm 


5 fm 


25 fm 


1 fm 


5 f m 


25 fm 




0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


q 


1.00 


28.7 


1.61 


2.80 


0.293 


0.249 


0.060 


0.702 


0.343 


0.393 


0.225 


0.216 


0.153 


54.5 


6.24 


9.46 


1.72 


1.52 


0.513 


1.10 


22.4 


1.78 


2.19 


0.285 


0.196 


0.053 


0.556 


0.296 


0.311 


0.189 


0.171 


0.124 


53.7 


8.00 


9.34 


2.01 


1.51 


0.553 


Initial condition: e' m '=28.7 GeV/fm fixed 




e q [GcV/fm ;i ] 


T q [GeV] 


s q [1/W] 


r 


1 fm 


5 fm 


25 fm 


1 fm 


5 fm 


25 fm 


1 fm 


5 f m 


25 fm 




0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


0.0 


3.0 


q 


1.00 


28.7 


1.61 


2.80 


0.293 


0.249 


0.060 


0.702 


0.343 


0.393 


0.225 


0.216 


0.153 


54.5 


6.24 


9.46 


1.72 


1.52 


0.513 


1.10 


28.7 


1.53 


2.72 


0.272 


0.235 


0.057 


0.591 


0.285 


0.329 


0.186 


0.179 


0.127 


64.7 


7.14 


11.0 


1.94 


1.73 


0.587 




n fl n 




FIG. 4: Profiles of the energy density e q (panels (a) and (b)), entropy density s q (panels (c) and (d)) and temperature T q 
(panels (e) and (f)) as functions of r\ calculated for different proper times r and different nonextensivity parameters q by using 
initial conditions with fixed e ( - m ' ) =28.7 GeV/fm 3 . We also obtain similar results for e q ,T q and s q with a — 1.25 fixed initial 
condition case. 



the fluid rapidity a q grows during the hydrodynamical usual hydrodynamic expansion (i.e., a q -\ > a q> \). 

expansion from its initial value a(rn,r?) = n. One can . ,,. , , , , 

, , . , , , , , lo summarize this part: one observes that nonexten- 

obscrvc that a lor q > 1 is decelerated compared to the a ., , a ., ^ , , , 

q H ^ sive huid (g-fiuid with q > 1) evolves more slowly than 
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the ideal fluid (with q = 1) 



where is the normal covariant vector of the isother- 
mals, 



E. Freezeout surface and single particle spectra 

We now present examples of single particle spectra 
emerging from our approach. We shall follow the sim- 
plest possibility in which they are expressed as an inte- 
gral of the phase space particle density over a freeze-out 
surface Ef |47j . 



E- 



^d 3 N 
dp 3 



d 3 N 



m^draTdydcj) 
9 



da^x) p^f e q(x,p). 



(27T)" JJk 

In the T-r) metric the surface element of Ef given by 



(33) 



dcTfi = (d<r T , da n ) — At rdrj I 1, 



(34) 



(n T ,n rl ) = 



dT 
Ifr 



dT 
drj 



(35) 



and At is the transverse area of the generated fluid. In all 
examples of applications to Au+Au collisions discussed 
in this paper we use At = 6.5 fm. The momentum of 
the produced particle in the r-r\ metric is given by 



77iT cosh(y — 77), —uit sinh(y — r\) 

T 



(36) 



where y is the observed rapidity (after freezeout). Using 
these expressions, the single particle density is given by: 



E- 



d 3 N 



d 3 N 



dp 3 rriTdmTdy 47T 



' ! | 'drj r f (?/) 



rriT cosh(y — if) ^|iiit sinh(y — rj) 

Tn T [rj) 



f q {y,v), 



(37) 



where 

f q {y,v) = 



l-(l-q) 



ttit cosh(y — rj) 



(38) 



and Tp is the freezeout temperature which is given by 
the corresponding freezeout energy density £f- In prin- 
ciple, the freezeout surface can be defined either as the 
surface of constant temperature Tp, or as the surface of 
constant energy density £f, or, finally, as the surface of 
constant entropy density sp (cf . Table |III|) . They all co- 
incide in the usual extensive case (9=1). In Fig. [S] we 
show as example freezeout surfaces (calculated for differ- 
ent values of parameter q and for different initial condi- 
tions) for Tp = 100 MeV. One observes a quite strong 
q dependence of the freezeout surface characteristics on 
these parameters. These dependencies are much weaker 
when calculated for surface of constant energy density 
and even weaker for constant entropy density (not shown 
here explicitly) . Note that values of Tp corresponding to 
freezeout conditions set by fixing £p or sf now depend 
on the parameter q (see Table UllJ). 

In Fig. [7] we show examples of single particle rapidity 
and transverse momentum spectra calculated for both 
types of initial conditions using Tp=100 MeV. Note that 
different types of the freezeout surface used are connected 
with using different sets of parameter (q, Tp), cf. Ta- 
ble lllip . Both distributions are sensitive to q, however, 
in the case of dN/dy this dependence is almost entirely 
due to the q dependence of the initial entropy density in 



0.6 , E (in) =28 .7 [GeV/fm*] » «n 






x=1.2fm 


- 9 = 1.10 


t =2.0fm 


■u-5.0fm 




t =25fm 




- - t =75fm 


1 \ 


■ "« 


» ! 


^....TImm^i 



s ? [GeV/fm 3 ] 



s q [GeV/fm 5 ] 



FIG. 5: (a) - evolution of the fluid rapidity a q (presented as 
a q — 77) as a function of the energy density e q for different 
values of r and for different nonextensivity parameters q with 
fixed e (in) = 28.7 GeV/fm 3 . (b) 
using a = 1.25 fixed initial condition case. 



the central region observed in Fig. [3] (panels (c) and (f)) 
and practically vanishes in the case of normalized distri- 
butions calculated for the constant £p freezeout surface, 
see the panels of Fig. [S] (a) and (b). This is because 
of the observed q dependence of the corresponding to- 
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TABLE III: The values of the freezeout temperatures Tf (in MeV) for different freezeout (F.O.) conditions used and different 
values of q investigated. 



q 


1.00 


1.05 


1.10 


1.15 


F.O.= T P fixing 


100 


100 


100 


100 


F.O.= £p fixing 


100 


91.8 


83.2 


74.3 


F.O.= sf fixing 


100 


89.3 


78.5 


67.4 




FIG. 6: Examples of freezeout surfaces of constant tem- 
perature Tf=100 MeV calculated for different values of pa- 
rameter q and for different initial conditions: with con- 
stant a = 1.25 (a) and with constant energy density 
£ F =9.82xl0- 3 GeV/fm 3 (b). 



tal multiplicities and is connected with the increase of 
the entropy observed in nonextensive processes, see the 
panel Fig. [5] (c) . We shall discuss this point in more de- 
tail below in Section [V] The weak residual q dependence 
observed in this case can be attributed to the (appar- 
ently very weak) effects of the EoS and freezeout surface. 
In what concerns px-spectra shown there for different 
initial conditions and freezeout surfaces, one observes a 
very strong dependence on q, which changes the slope of 
Pt considerably. It is interesting to note that, as seen 
in Fig. [3 the px distributions apparently are not sensi- 
tive neither to the type of initial conditions nor to the 
freezeout surfaces used. 



In the pt distributions, the slope depends on both q 
and Tp and increasing Tp while keeping constant q gives 
a similar effect as increasing q at fixed Tp. On the whole 
one observes tendency that transverse expansion as mea- 
sured by these distributions gets stronger with increasing 
nonextensivity, i.e., with increasing q. 




FIG. 7: dN/dy and pt spectra (panels (a)-(b) and (c)-(d), 
respectively) obtained from the g-hydrodynamical evolution 
with constant Tf =100 MeV and for different values of pa- 
rameter q and for both types of initial conditions: with fixed 
a = 1.25 (panels (a) and (c)) and for fixed e (in) = 27.8 
GeV/fm 3 (panels (b) and (d) ). Rapidity spectra are obtained 
by integrating Eq. (|37[) over pr £ (0,6.0) GeV/c whereas pt 
spectra are obtained by integrating Eq. (|37|l over \y\ < 0.5. 



IV. COMPARISON WITH EXPERIMENTAL 
DATA 

We shall confront now our approach with experimental 
data. Because of still explanatory character of our work 
we limit ourselves only to comparison with some selected 
rapidity and px distributions. At this stage no attempts 
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£ (in) =28 ? [Qgv/fm*] 



s F =9.82xlO~ 3 GeV/fm 3 surface 



1.00 



1.05 
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1.15 



FIG. 8: (a) and (b) - normalized rapidity distributions 
\/N dN/dy (defined in the same way is in Fig. 0) calculated 
for constant energy density freezeout surface £f=9.82 -10~ 3 
GeV/fm 3 for different values of q using initial conditions with, 
respectively, fixed a = 1.25 and fixed e (in) = 28.7 GeV/fm 3 . 
Notice that there is only very weak q dependence confined to 
small and large regions of rapidity y. (c) - q dependence of the 
total multiplicity V to t obtained from the g-hydrodynamical 
evolution with fixed e'" 1 ' = 28.7 GeV/fm 3 initial condition 
and e F = 9.82 ■ 10~ 3 GeV/fm 3 freezeout condition. The total 
multiplicity TVtot increases linearly with q. 



for exact fits have been made. They must wait for a 
more detailed version which, for example, would account 
for the possible changes of the nonextensivity parame- 
ter q during the collision process as mentioned in Section 
HI The same remarks apply to the potentially promising 
analysis of anisotropic flow or particle interferometry (for 
example, in the way as it was done in [33l l34l. l35l] ) . which 
we postpone until 1 + 2 dimensional version of our ap- 
proach accounting for expansion in transverse directions 
will be available in the future). Because, as was shown 
in Section fill El the most sensitive for q dependence are 
Pt distributions, we start with them and show in Fig. [9] 
that data by Q prefer q = 1.08 and Tp = 100 MeV (we 



attribute the visible discrepancy at largest values of pr to 
the contamination from quark jets which carry large mo- 
mentum in the initial stage of nuclear collisions and which 
are not accounted for g-hydrodynamical model. With 
these values of q and Tf data provided by [49j for dN/dy 
distributions and by [48| for pr where compared with 
predictions of different initial conditions characterized by 
e' ln ) , see Fig. [TU] As one can see, the g-hydrodynamical 
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FIG. 9: Comparison of q-hydrodynamic model with experi- 
mental data observed by STAR collaboration 48] performed 
using a = 1.30 and T F = 100 MeV for q = 1.08, 1.09 and 
1.1 values (with corresponding values of e' m ' = 21.2, 205 and 
19.7 GeV/fm 3 ). The best agreement is obtained for q — 1.08. 



model with g-Gaussian initial condition can reproduce 
reasonably well both the rapidity and transverse mo- 
mentum distribution data simultaneously. It should be 
stressed at this point that with the parameter q > 1, 
which according to the general philosophy of the nonex- 
tensive approach accounts for all possible intrinsic fluctu- 
ations in the system @, H, Q , our model also accounts for 
the possible presence of resonances [H, EH which there- 
fore, to avoid double counting, should not be added inde- 
pendently. It must be noticed that in the present version 
we do not, in fact, account for the possible creation of a 
quark-gluon plasma (QGP) phase. For this one should 
use a more elaborate version of EoS than discussed here 
in Section Till Bl Notwithstanding all that, one can say 
that a simple q-hydrodynamical model reproduces exper- 
imental data reasonably well using e^ 111 ) = 19.0 -H 22.3 
GeV/fm 3 (cr = 1.28 1.32), T F = 100 ~ 120 MeV, and 
q = 1.07 -=-1.08. 
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FIG. 10: Comparison of the g-hydrodynamic model with experimental data on rapidity [40] and pt [481 ] distributions calculated 
for g=1.08 (fixed) and Tf=100 MeV, as in Fig. [9l but for different (Gaussian) initial conditions parametrized by e' m '. 



DISCUSSION: CAN PERFECT 
(/-HYDRODYNAMICS MIMIC 
d-HYDRO DYNAMICS? 



Nonextensive / dissipative correspondence 
formulation 



Our starting point is observation made at the end of 
Section IIII Dl that g-fluid evolves more slowly than an 
ideal fluid. To this one can add observation from Sec- 
tion IIII El above that transverse expansion measured by 
the behavior of px-spectra is much stronger in g-fluid. 
Those are precisely features observed in viscous fluids 
(cf., for example, [35|]). Let us then treat these obser- 
vation seriously and look more closely for the possible 
connections between g-fluid and viscous fluid apparently 
emerging from our g-hydrodynamical model. 

Let us start with reminding the possible physical mean- 
ing of perfect g-hydrodynamics. It originated from the 
modified Boltzmann kinetic equation (|3a|) in which a new, 
(/-generalized, version of the Boltzmann molecular chaos 
hypothesis @, 0, HI, H, E3, H3| has been used in the form 
of Eq. Q. It can be introduced in different ways [66| 
but effectively it always amounts to postulating a new 
kind of equilibrium, which includes some interactions and 
in which some stationary state is formed (8j summarily 
characterized by parameter q. In our case it leads to 
Eq. ([20]) . which is formally identical to perfect hydrody- 
namical equation but with all usual ingredients replaced 
by their g-counterparts (perfect means here that there is 
nothing on the r.h.s. of Eq. (120p ). It is natural to ask 



how Eq. (|20|) would look like when written in terms of 
the usual perfect hydrodynamic (with q = 1) and some 
reminder depending on the parameter q. Because, as we 
have seen, in general q differs only slightly from unity, 
q — 1 <C 1, it is tempting to simply expand Eq. (|2"0|) in 
small parameter \q — 1| [ill [Hoj . However, as shown in 
Appendix[C] in such case one faces some unsurmountable 
problems because terms multiplying \q— 1| are not small 
enough in the whole of phase space. We shall therefore 
follow more general approach. 

All our results presented above come from the Eq. 
(I20p , which is equation for perfect g-hydrodynamics. No- 
tice that nonextensivity affects not only the thermody- 
namical quantities like energy density e and pressure P 
but also the flow velocity field u M (a;): 



e(T) 
P(T) 
u»{x) 



e q (T q ) =e(T q ) + As q (T q ), 
P q (T g ) = P(T q )+AP q (T g ), 
u£0) = u"(a:) + Su^{x). 



where (x) is formally a solution of the equation which 
has form of dissipative hydrodynamical equationJH, [24], 

HE HI HE Hi Ha H M Hi , 



;W - PA"" + 2W^u v) + tt^ 



(39) 
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The notation used is: 



3n, p = p q + n, 



W^ = w q [l + 1 ) A^u* 
W^W V 



W q [l+~f} 2 

w q SufSuf 



(40a) 
(40b) 



(40c) 



where i is energy density, P pressure, energy or heat 
flow vector, tt^ the (symmetric and traceless) shear pres- 
sure tensor and where 



w = E 



Pa 



7 = u^Su^ = --6u qli 5u%, 



(41) 
(42) 



and 

a<^b v> 
whereas 



[~(A^ + A£A£)-~A^A ACT ]aV 



n 



r u> 9 [ 7 2 + 2 7 ]. 



(43) 



This last quantity can be regarded as a bulk pressure to 
be used below. 

Now comes crucial point of our argumentation. To 
proceed further we shall assume that there exists some 
temperature T and velocity field Su^ satisfying the fol- 
lowing relations: 



P(T) = P q (T q ) 
e{T) = e q {T q ) + m 



(44a) 
(44b) 



(e and P are energy density and pressure defined in the 
usual Boltzmann-Gibbs statistics, i.e., for q = 1). In this 
case one can transform equation (|39[) into the following 
equation, 



e(T)u»u u -{P{T) + n)A' ll/ +2W^ 



= 0,(45) 



which has the familiar form of the usual d- 
hydrodynamical equation. This means that perfect 
g-hydrodynamics represented by Eq. (|20j) can be re- 
garded as being formally equivalent to some form of 
d- hydrodynamics as represented by Eq. (|45|) . We shall 
call this observation the NexDC correspondence (and, 
respectively, we shall call Eq. (|4*4")l with Eq. (|40bp and 
(|40c[) the NexDC relations). This observation can be 
traced back to the fact of generic non-conservation of 
global entropy in nonextensive systems, cf. Eq. (|2ip . 
visualized in Fig. [5] as increase of the multiplicity with 
increasing q. 



B. Nonextensive/dissipative correspondence - 
consequences 

We shall now present shortly the most specific immedi- 
ate consequences of NexDC correspondence: the entropy 
production and estimations of the corresponding trans- 
port coefficients. 



1. Entropy production in q-hydrodynamics 

Let us start with observation that Eq. (|43[) and NexDC 
relations (|44|) lead following form of gr-enthalpy, 



e q (T q ) + P q (T q ) 



e{T) + P(T) 
[1 + 7] 2 : 



(46) 



which can be also used in definition of 7 because w = 
Ts = e + P and 1/(7 + 1) = y/l - 3U/w (s is the entropy 
density in the usual Boltzmann-Gibbs statistics). Notice 
that in NexDC one has that 



WW f ,= -3ILw, 

■k» v w v = -211 w* 4 , 
7r u ^=6n 2 . 



(47a) 
(47b) 
(47c) 



Suppose now that we define a true equilibrium state as a 
state with q = 1, i.e., with no residual correlations be- 
tween fluid elements and no intrinsic fluctuations present, 
with energy momentum tensor 



T£ = T" u = e{T)u^u v - P(T)A^ 



(48) 



and with equilibrium distribution given by the usual 
Boltzmann distribution, 



f eq (x,p) = exp 



p^u^x) 



kB,T(x) 



(49) 



In this case the state characterized by f q (x,p) given 
by Eq. (TT4"]) must be regarded as some stationary state 
existing near equilibrium. Therefore, because we ex- 
pect that \q — 1| is small, we can define a near equi- 
librium state defined by the correlation function h q in 
Eq. (0J for which the energy momentum tensor is T q ^ u = 
{e q + P q )u^u v q - P q g^ v , cf., Eq. It means then that 

we can write 



where 



ST»" = -UA^ + W^u" + W v u ii + n^. 
Using now Eq. (|4~4"|) we obtain the relation 



where 



5c 



j=y/l + 5e q -l, 

_ e(T) - e q (T q ) 
1 e q (T q ) + P q (T q y 



(50) 
(51) 
(52a) 

(52b) 
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which connects the velocity field u (solution of the dis- 
sipative hydrodynamics given by Eq. (|45[) 1 with the ve- 
locity field u q (solution of the q-hydrodynamics given by 
Eq. (HOD). 

(1 + 1) dimensional case discussed here, 
can always parameterize these velocity fields 



In 

one 

by using the respective fluid rapidities a q and a, 

u%(x) 

[cosh (a — rj) , \ sinh (a — n)~\ . 
cosh(a g — a) — 1 one has that 



[cosh (a q — rf) , ^ sinh (a q — ??)] and u^{x) 
Because 7 



cosh(a g — a) = y/l + 5e 



1 1 



(53) 



which provides us with a connection between u and u q . 
From Eq. (|53p. one obtains finally 



log (e q + yr" 



Se, 



(54) 



We abandon here another solution of Eq. (|53|) . namely 
that a — a q + log (e q + + Se q j , because it leads to 
the entropy reduction, i.e., for it [su M ] ;M < 0, for q > 1. 
Taking the covariant derivative of Eq. (|50[) and multiply- 
ing it by u v we obtain 



0. 



(55) 



Therefore, although in ideal g-hydrodynamics the q- 
entropy is conserved, i.e., [s q u^] [f j_ = 0, we can rewrite 
it in the form corresponding to dissipative fluid with en- 
tropy production, 
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(56) 



To illustrate this we show in Fig.[TT]the expected entropy 
production as given by Eq. (|56[) . Notice that su^ > 
for the large 77 region at any r (but especially for the 
early stage of the hydrodynamical evolution). It supports 
therefore a dissipative analogy of the g-hydrodynamics 
mentioned before and leads us to very interesting con- 
clusion that equilibrium state which is generated in the 
high-energy heavy-ion collisions may in fact be the q- 
equilibrium state which can be regarded as some station- 
ary state near the usual (i.e., q=l) equilibrium state and 
which contains also some dissipative phenomena. 



2. Calculation of transport coefficients from 
q-hydrodynamics 

There are different formulations of d-hydrodynamics 
0, Hi, H, El, [13, H, H, HU HI, HI . In what follows we 
shall choose for comparison only the 2nd order theory of 
dissipative fluids (in particular as given by [2J, |25|) leav- 
ing investigations of other approaches from nonextensive 
perspectives for future investigations. As is known, this 
theory does not violate causality (at least not violate the 
global causality over distant scale given by the relaxation 
time) , on the other hand it contains now some dissipative 
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FIG. 11: Evolution of the entropy production [su^]-^ (panels 
(a)-(c)) and ratio [su M ] ;(J /s (panels (d)-(f)) as function of n for 
different values of r for q = 1.08 with e (in) = 22.3 GeV/fm 3 
(or a = 1.28). The error bar is estimated by the value of 
\sqUg]^ obtained in the numerical calculation which should 



be zero in an analytical calculation. 



fluxes like heat conductivity, bulk and shear viscosity. We 
shall now see to what extent these transport coefficients 
can be calculated in g-hydrodynamics. 

To this end let us start from considering more closely 
respective entropies. Dissipation is connected with the 
production of entropy and in [24], [25[ the most general 
off-equilibrium four-entropy current er^ can be written 
as 



P(T)f3» + v (TJi v + ST"") + Q 



(57) 



where fi» = /T and Q» = (ST^) is some function 
which characterizes the off-equilibrium state. In the case 
of the g-entropy current (0 the NexDC conjecture (i.e., 
Eqs. (|40b[) and (|4"6")) ) leads to the following off-equilibrium 
state: 



Q M = Q» = 



T 



(58) 
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(with x = t; ]/ 1 " 



311 



1) which results in 



+ — +x|s« M + — y (59) 

Notice that, because of the strict g-entropy conservation 
assumed here, using Q M = one always gets = 
0. It means that, although there is no production of q- 
cntropy, there is production of the usual entropy, i.e., our 
g-system is really dissipative in the usual meaning of this 
word. 

Let us now be more specific and use the most general 
algebraic form of Q M , calculated up to the second order 
in the dissipative flux, as given by (25| 



?2nd 



2T 



(60) 



T T 

Here /3i=i,2,3 are the corresponding thermodynamic co- 
efficients for the, respectively, scalar, vector and tensor 
dissipative contributions to the entropy current whereas 
ai = o,i are the corresponding viscous/heat coupling co- 
efficients. The n is bulk pressure defined before in Eq. 
[67j . In the NexDC one has that 

r 2n d SU^ + Ti s t- 



?2nd 



T 



where 



2nd 



2w 



T lst = -(a + 2ai)n. 



(61) 

(62a) 
(62b) 



Q M can be then expressed by polynomials in the bulk 
pressure II defined by Eq. (|43)l . It is then natural to ex- 
pect that the most general entropy current in the NexDC 
approach has form: 



# uU = r(n)^ + T(n)^, 



(63) 



where T, T are (in general infinite) series in powers of the 
bulk pressure II. In this sense the can be regarded 
as the full order dissipative current. 

In general one has entropy production/reduction, 
att 7^ 0, however in case when r(n) = T(II) = x one 
has = so one can write the full order dissipative 
entropy current as 



'full 



( X + £)su>* + {x - 



where T and T are determined by x — (r - 
£ = (F - T)/2. From two solution for (r, T), 



(64) 
T)/2 and 



T T ( 311 \ T T-T a , . 

5-e V 1 - — 1 ■ 2 s ir^ a 



only (|65ap is acceptable because only for it w u (3^ lU < 
(i.e., entropy is maximal in the equilibrium [25j and this 
is because (T — T q )/T q is always positive for q > 1). 
In this way we finally arrive at the following possible 
expression for the full order dissipative entropy current 
in the NexDC approach: 



'full 



2T 



1 - 



1-™ 

w 



2(T-T q ) W 1 



T 



(66) 



Limiting ourselves to situations when T/T q w 1 and ne- 
glecting terms higher than 0(3II/-u;) 2 , one obtains that 



Q* 1 



full 













\ w ) 



(67) 



Comparing now Eqs. (|rJ2^1 and (|6T|) one gets that (66 



/?! = -, 



A)- 



6/?2 = 5-, 
2iu 



a + la x = 0. (68) 



Since in the Israel-Stewart theory [24| the relaxation time 
r is proportional to thermodynamical coefficients /?o,i,2> 
it is naturally to assume that in our NexDC case r oc 1 /w, 
i.e., it is proportional to the inverse of the enthalpy (no- 
tice that for classical Boltzmann gas of massless particles 
one obtains & = 3/w (25ll3^). 

We shall now derive the bulk and shear viscosities 
emerging from the NexDC approach. Let us start with 
observation that the local entropy production by the full 
order entropy current Eq. (|66[) can be also written as 



full j 



where $ M = 



= [(1 + X)*"]i„ 



and ^ = su^ 



T 



(69) 



Because 



conservation of g-entropy, = 0, is equivalent to [(1 - 
x)^*];m = 0' therefore using Eq. (T4"6l one gets that 



W V W V 
3ILT 



2UT 1 



(70) 



and (see Appendix [Dlfor details of derivation of Eqs. (fTTj) 
and 1(731 ) 



'full;/j ' 



-.--(wu"Xj-—Y lt + —Z lu> , 



where 



n 



T 



^ui/ — 77 



5^ 



(71) 

(72a) 
(72b) 
(72c) 



and 



2n ' 



(72d) 
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One can now use Eq. (|47|) to eliminate term proportional 



to heat flow 



In this way one avoids explicit contribu- 



tion to entropy production from the heat flow ^L- , which 
is present in Eq. (|71[) when one discuss baryon free fluid, 
in which case the necessity to use Landau frame would 
appear. As one can see Eq. (|69[) is covariant and there- 
fore it does not depend on the frame used. After that 
one obtains that 



'full;/; 



(73) 



Eq. dZ2 

viscosities from 



can be now used to find the bulk and shear 
given by Eq. ([6"5)) . The positive 
transport coefficients, bulk viscosity C and shear viscosity 
77, can be estimated by writing the entropy production 
Ofjjj.u in the following form: 



'full;;/ 



n 2 



1 fj,is 



2r]T 



> 



(74) 



and using Eq. (|47[) . We arrive then at the sum rule con- 
necting transport coefficients (expressed as ratios of bulk 
and shear viscosities over the entropy density s), 



1 _3_ 

C/s 77/* 



'full; 



n 2 



(75) 



This is as far as we can go. The heat conductivity, as 
shown above, can be expressed by two other transport 
coefficients for which we have only one equation in the 
form of sum rule (|75|) . To proceed any further and to 
disentangle ([75)1 one has to add some additional input. 
Suppose then that we are interested in the extremal situ- 
ation, when total entropy is generated by action of shear 
viscosity only. In this case one can rewrite Eq. ([75)) as 



resulting in 



7T^ 



611 



(wu x X x ) 



v 7(7 + 2) 



s 3(7 + 1) 2 



n t 



SU X Xy 



(76) 



(77) 



Note that Eq. ([77)) allows all values of 77/s, in partic- 
ular that 77/s < j-, what violates the limit obtained 
from AdS/CFT correspondence that 77/s > 1/Att [36| . 
To impose this limit we shall now use Eq. ([75]) . This can 
be done only in the region where r.h.s. of Eq. ([77)) is 
smaller than (or equal to) 1/Att (in which case we put 
77/s = 1/4-71"), otherwise, because of our earlier assumed 
limitation we put Q/ s = and use Eq. flTTJ to evalu- 
ate 77/ s. The corresponding results for £/s and 77/s are 
shown in Fig[T2] (a) and (b), respectively. Notice that 
when the r.h.s of Eq. ([77)) approaches 1/47T, (/s given 
by Eq. (j75[) approaches infinity. All curves presented in 
Fig. [T2]were calculated for Au+Au collisions with iden- 
tical set of parameter as the best fit presented in Fig. ITD1 
above. 
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FIG. 12: The NexDC predictions for the ratios of bulk (panel 
(a)) and shear (panel (b)) viscosities over the entropy density, 
C/s and rj/s, as function of temperature T and calculated for 
a number of space-time rapidities y = | In j±r (the same as 
in Eq. (|24b[) using g-hydrodynamical model). 



VI. SUMMARY AND CONCLUSIONS 

We have presented a nonextensive version of the hy- 
drodynamical model for multiparticle production pro- 
cesses, the q-hydrodynamical model, which is based on 
the nonextensive statistics represented by Tsallis entropy 
and indexed by the nonextensivity parameter q. In do- 
ing so we have followed the usual approach originating in 
the appropriate kinetic equations formulated in nonex- 
tensive form in (l7| . We have found the nonextensive 
entropy current which satisfies not only the nonextensive 
if-theorem, Eq. ([5]), but also the q- version of thermo- 
dynamical relations Eq. (|L7[) . The (1+1) dimensional 
g-hydrodynamics with the g-Gaussian initial condition 
and the g-EoS can reasonably reproduce the single parti- 
cle spectra observed at RHIC energy for q = 1.07 -j- 1.08 
for T F = 100 120 MeV if quark jet contributions to 
Pt spectra are small, i.e., up to transverse momentum 
range around pr < 6.0 GeV/c. We also found a pos- 
sible correspondence between the g-hydrodynamics and 
the usual (q = 1) d- hydrodynamics (NexDC) as provided 
by the Eq. (gSJ) with NexDC relations Eq. (JUJ). Based 
on this correspondence, we have evaluated the entropy 
production in relativistic heavy-ion collisions at RHIC 
energy using results of our perfect ^-hydrodynamics (un- 
derstood as an approach without any g-viscosity effects 
added). The fact that when comparing with data one 
finds that q > 1 means that one, indeed, has some dy- 
namic factors present, detailed form of which is not yet 
disclosed but which summarily can be accounted for by 
the nonextensive approach and which action is summa- 
rized by the parameter q — 1 . 

In what concerns the obtained pt dependence, our for- 
mula continues attempts made to interpret power-law 
spectra as new kind of equilibrium phenomena for the 
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whole pt range pushing the usual interpretation via on- 
set of "hard" collisions (imposed on the "soft" ones) to 
really high values of p T (cf., |H[l| and references 
therein). In such approach there is no characteristic 
scale at which the transition from " soft" (or locally ther- 
malized) to "hard" (or unthermalized) dynamics occurs 
which appears in the conventional descriptions using vis- 
cous hydrodynamics, as for example in [34], |51( . 

One of the results of our investigation is that fluc- 
tuations in the initial conditions seem to be the most 
important part of the hydrodynamical model, which by 
using Tsallis statistics, attempts to account for any pos- 
sible fluctuations in some general, model independent, 
way. This is quite reasonable result because at the ini- 
tial stage our system consists of a relatively small num- 
ber of degrees of freedom and is therefore more sensitive 
to any fluctuations. On the contrary, at freezeout this 
number is much bigger and the system is only weakly re- 
sponding to any fluctuations. This finding agrees nicely 
with recent analysis of the elliptical flow performed by 
using a hydrodynamical approach in which an attempt 
was made to account for fluctuations (without, however, 
using q-statistics), see [52[. On the other hand, however, 
it should be remembered that analysis presented here is 
considerably simplified by using the same nonextensivity 
parameter q at all stages of the collision process. There 
is therefore room for improvements which will facilitate 
comparison with data. One can argue that intrinsic fluc- 
tuations existing in different stages of the collision pro- 
cess are of different (albeit connected) dynamical origin 
and therefore parameters q for the initial conditions, for 
the EoS and, finally, for the hydrodynamical expansion 
should be allowed to have different values (and should 
be also different for the longitudinal and transverse dy- 
namics). The other problem would be how to connect 
our q parameter expressing fluctuations with fluctuations 
in all momentum observables as seen when analyzing a 
non-ideal liquid as it was done, for example, in [53j |. We 
plan to address these problems elsewhere. In any case, 
similarly to the fact that concept of ideal fluid is never 
realized in nature [3a ] (the bound of rj/s > 1/4% found 
in [36| being strong argument supporting this), the q = 1 
case should be replaced by investigations of the (/-fluid 
with q > 1. 

In this context the natural question arises concern- 
ing the deeper physical meaning of the q-hydrodynamic 
proposed here. The most important observation dis- 



cussed in Section[V]is the apparent correspondence found 
between the perfect (/-hydrodynamics and the usual d- 
hydrodynamics, which we call the NexDC correspon- 
dence. It allows calculation of transport coefficients of 
viscous fluid in terms of parameters of q- (ideal) fluid, i.e., 
essentially as dependent on single parameter which, as it 
was already stressed many times, represents summary ef- 
fect of many possible dynamical factors, without entering 
into dynamical details (i.e., in a purely phenomenological 
way). The detail discussion of the NexDC phenomenon 
is, however, outside the scope of present paper and we 
plan to address is elsewhere. 

We close by remark that hydrodynamics can also 
be derived using Information Theory with its method 
of maximization of information entropy under some 
specific constraints 54]. It is therefore plausible that 
our results could also be derived using a nonextensive 
version of Information Theory (in the way as [l3l . |55| 
can be regarded as a nonextensive generalization of 
the Information Theory approach to single particle 
distributions obtained in the multiparticle production 
processes proposed in [HI). We shall not pursue this 
possibility here. 
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APPENDIX A: DERIVATION OF EQS.([26]) AND 

Consider some details of the (1 + 1) dimensional rcl- 
ativistic hydrodynamics under assumption that one ig- 
nores the transverse expansion of the fluid. With — 
diag(l, — t ) and four fluid velocity — [cosh(a g — 
77), — Tsinh(a 9 — 77)] the projection matrix is equal to: 



— sinh 2 (a g — 77) 
t cosh(a g — 77) sinh(a ? — if) 



t cosh(a 9 — 77) sinh(a (? — 77) 
-t 2 cosh 2 (a, — 77) 



(Al) 



The nonvanishing components of Christoffel symbols are derivative of fluid velocity, which is defined by u^ v — 
= TJ' = 1/r and — r, therefore covariant 
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<9„u M + r^u A , has the following form: 

, da„ 



da n 



u q . T = sinh(a 9 - r >)~Q^-> u q\r = ~ cosh(a 9 - ^ g > 
u T q . v = sinh(a g - 77) ^p, vP q . n = ~ cosh(a g - r?)^. 
Using these expressions, one obtains 



cosh(a g — 77) sinh(a g — rj)—^- 

OT 

1 ■ 1 2 / \ da„ 

I — smh (a q - ri)—-, 

T Of] 

-rcosh 2 (a„ - r\)%r^- 

OT 

— cosh(a„ — 77) sinh(a„ — 77)— -2-. 

OTj 

Because g?" = for /1, v = r and 77, Eq. (f2"3")l is reduced 
to the following two equations: 



?7 M A V V 



dP 

-A —± 
qTT dr 



1 dP 

-A qTV ^=0, (A2) 



{s q +P q ){u> i q & mv u u q . ll } 

A dP q 1 A 

T 



$77 



5P„ 



OTT 5r ^ r 2^™ 577 



0. (A3) 



The equations (|A2p and (|A3j) are equivalent, therefore 
one has only one equation, 



tanh(a g — 77) cta^ 
r 9?7 
1 <9P, 



<9P , ,,, 

. tanh(a ,_ f? )_a + __^ = 0> (A4) 
I 



which is the Eq. (|27)) . Since the four divergence of the 
fluid velocity is given by 



- smh a, - 77)— -2- + - cosh(a, z - 77)— — , 
or t or/ 



the Eq. ([23)1 can be written as 



<9e g sinh(a„ - 77) de q 

cosh(a - 77)— - H 2 ^ 

or r a?7 

f 

+ (e g + P q )i sinh(a 9 - rj)-^ 

+ i cosh(a 9 - 77) ^ I = 0, 



leading immediately to 



de q tanh(a g — 77) de q 
dr t 



drj 



. f , , , da„ 1 da a 1 

+ (e, + P 9 ) |tanh(a 9 - 77)^ + = 0, 



(A5) 



which is Eq. ([26 



APPENDIX B: NUMERICAL METHOD USED 



For the purpose of numerical calculations we express Eqs. 

(n+i) 1 r (™) l (™) 



and (|27p in the form of finite difference equations 



(n) _ (n) 

(n) I £ g(j+l) £ g(j-l) 




= 



(Bl) 



and 




(B2) 
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The subscript (J) and superscript (n) represent the corresponding grid number in the r? and r space with grid spacings 
Ar? and At, respectively, i.e., with rjj = j'Ar? and r„ = r + nAr. The coefficients appearing in the above equations 
are defined in the following way: 
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where c s is function of eJ") (due to the equation of state P q (e q )), one can rewrite above two equations in the 
following form: 
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Eliminating a 
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from the above two equations, one obtains 
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where relations a 
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APPENDIX C: INADEQUACY OF THE SIMPLE 
EXPANSION IN |g-l| 

From the previous experience in appl ying ^statistics 
to multiparticle production [13, HH, [S Elt 0, EH we 
know that \q — 1| < 1. It seems then natural to argue 
that (see, for example, [ill. [Hoi]) one could simply expand 

1/(1-9) 



f q (x,p) = l-(l-q) t 



<T(x) 



from Eq. (T 



z = 1 — q, retaining only terms linear in z and get: 

1 



f q (x,p) =f(z) = [1-z- A\*i 



■ ln(l - z ■ A) 



f(z = 0) + z 



df(z) 



dz 



(CI) 



(here A = A(x,p) = p and the arguments (x,p) are 

suppressed for clarity). 

However, such expansion can only be performed under 
some conditions, which we shall clarify in what follows. 
It is straightforward to show that to gets first step of 
expansion, 



exp 



-A ■ 



A 



exp [— A] ■ exp 



.4 2 



it is necessary that 



z ■ A(x,p) < 1. 



(C2) 



(C3) 



The second step needed is to additionally expand the 
exponent and this requires 



z ■ A 2 (x,p) < 2. 



(C4) 



When this is satisfied, one finally gets f q (x,p) in terms 
of f q=1 (x,p) only, 



fq(x,p) ~ fq=l(x,p) + (l-q) 



A 2 (x, P ) 



fq=l(x,p). 

(C5) 



At first this procedure looks very promising because 
using it one gets 



(C6) 



where T q —i is the usual energy-momentum tensor for the 
equilibrium of the Boltzmann-Gibbs statistics, i.e. the 
one usually used when describing an ideal fluid, 



~ (27 



— J— p»p exp(-_ ) 



(C7) 



whereas the nonextensive correction tensor t£ v is given 
by 



fj 



d 3 P 



uv _ 

q - {2ttYJ p° 
•exp(-— ) 



p»p v ■ 



p ■ u 
~~T~ 



1 /p ■ u^ 2 

2 V~ 



,(C8) 



However, in our case condition (|C4[) would impose too se- 
vere constraints on the allowed q and the region of phase 
space, p and x, considered rendering this approximation 
rather unpractical for our purposes. 



APPENDIX D: DERIVATION OF EQ. (71]) AND 
r)/s, EQ. (73]). 



The entropy production is given by 



where 



Then, 



<n- lf , = [^%, = ^ 1 ^ + ^ ( m ) 



3nr 2nr 

= ^ {2W - 3tV} . (D2) 
olli 



ft=f^ + ^] ^ + —ffl., (D3) 
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where t/^ = 2W v W v u» - . The tp^ is explicitly 

written as 

v/; = {2w v w u u" - sw^"}.^ 

= [AW^u* + 2W v u^ - 3< M ] W + [-3W„ ;/ J 

and 

-1 



6nT r,n 



-1 

TT 



l. 



i 



n 



3 

~ 7 F' 



(D4) 



Hence we obtain that 



n 



i 

n 
n 



T 



2 f;p 



9 u n a,T 



n 



7 * M ~ ^ T 



Finally one arrives at 



'fillip 



n I £ 



n t 

w 

n - w 



= -X„(«m")=-y„^ r + Z^^. (D6) 



This is Eq. <J7T]). Using now tt^W,, = -2nH^ one gets 
that 



'full;/* 



, „ N n - IE" tt^ 



2ht + ^ 



r 



(D7) 



(D5) which is Eq. ([75]). 
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We mention at this point that the so called kinematic 
thermalization used here in which equilibration of ener- 
gies is due to the collisions, has been recently contrasted 
with the so called stochastic thermalization based on the 
process of erasing of memory of the initial state resulting 
in a state of maximal entropy and coinciding with the 
above thermal equilibrium state, see [52] and references 
therein. 

[65] Notice that, in principle, both the jVp ar t and the AE 
are also fluctuating quantities but we shall not consider 
these fluctuations here. One can argue that they are to 
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initial conditions considered here. 
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analysis of the g-Hagedorn model in [501 ] . 

[67] Notice that whereas the time evolution of II is controlled 
by Q-hydrodynamics (via the respective time dependen- 
cies of e q , P q and x) its form is determined by the as- 
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tropy production in the standard 2 nd order hydrody- 
namic theory [24], [25| is never negative. 

[68] The fact that we obtained non-zero coefficients /3;=i,2,3 
and couplings a?;=o,i for dissipative flux of Eq. (|60p as 
found in Eq. (|68[) means that d- hydro dynamics obtained 
via NexDC from q-hydrodynamics accounts for all second 
order terms. One may conclude than it seems that such d- 
hydrodynamics with full order entropy current has global 
causality (over a distance scale given by relaxation time). 
However, question whether NexDC violates causality re- 
mains so far unsettled. 



